Imaging-based Parametric Resonance in an Optical Dipole Atom Trap 
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We report sensitive detection of parametric resonances in a high-density sample of ultracold Rb 
atoms confined to a far-off-resonance optical dipole trap. Fluorescence imaging of the expanded 
ultracold atom cloud after a period of parametric excitation shows significant modification of the 
atomic spatial distribution and has high sensitivity compared with traditional measurements of 
parametrically-driven trap loss. Using this approach, a significant shift of the parametric resonance 
frequency is observed, and attributed to the anharmonic shape of the dipole trap potential. 

PACS numbers; 37.10.Gh, 37.10.Jk 



Parametrically driven processes are ubiquitous in na- 
ture, and occur in a range of classical and quantum sys- 
tems [l|, 0, 0, 0- In a parametrically driven system, 
a system parameter may be harmonically varied in time. 
For a one-dimensional harmonic oscillator with a char- 
acteristic frequency ujo, such excitation generates reso- 
nances at frequencies usp = where n is an integer n 
= 1,2,.. A familiar example is the inverted pendulum, 
which may be stabilized against decay of small oscilla- 
tions by appropriately driving the pivot point. One im- 
portant role of parametric excitation in atomic physics 
is its influence on ultracold atoms confined in magneto- 
0£tical 0, ma gnet ic 0, S, S , and optical dipole traps 
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[ill . Il2l . Il3l . [ij, ll5| . For situations where long lived traps 



are desirable, excitation due to noise in the trap param- 
eters normally heats the atoms and thus leads to atoms 
being expelled from the trap. Under some circumstances 
parametric driving of a system with an anharmonic po- 
tential can also cool the atom cloud 0, H, @| ■ The lifetime 
of the trap is limited by technical noise in any of the trap 
13l . Il4l | , two important ones for optical 
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parameters 

dipole traps being the pointing stability and the intensity 
noise characteristics of the trapping lasers. In addition 
to this deleterious effect, parametric resonance is a useful 
and widely employed tool for characterizing the shape of 
the trap, including the harmonic frequencies of the trap 
[l^ . This is important in experiments where the spatial 
shape of the trapping potential is needed for interpreta- 
tion of the results. In experiments in our laboratory, we 
are interested in obtaining a high density of atoms in an 
optical dipole trap, and in knowing the peak density as 
well as possible. This requires that the number of atoms, 
atom temperature, and shape of the trapping potential 
be well known. A main goal of these experiments is to 
study the electromagnetic analog of Anderson Localiza- 
tion [H, [13, \m, 19] of hght in an ultracold atomic gas. 

In this paper we present experimental results on para- 
metric excitation of a sample of ultracold ^^i?6 atoms 
confined in an optical dipole trap. Such studies are fre- 
quently done by measuring the relative number of atoms 
that are removed from the trap as excitation parame- 
ters are varied. In the present case we measure instead 



modification of images of the ultracold cloud following 
parametric excitation and a subsequent period of free ex- 
pansion. Analysis of the images provides a more sensitive 
measure of the effects of parametric excitation, including 
clear evidence of the frequency shift of the parametric 
resonance frequency due to trap anharmonicity ^, tOi, |2Q[ . 
In the following we provide a brief description of the ex- 
perimental apparatus. This is followed by presentation 
of the measured images and their analysis. 
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FIG. 1: Schematic drawing of the experimental apparatus. 
In the figure QUEST stands for quasi electrostatic trap, and 
CCD represents charge coupled device. Drawing not to scale. 

A schematic of the experimental apparatus is shown in 
Fig. 1. In the figure, the central part of the experimen- 
tal apparatus is a magneto optical trap (MOT) confin- 
ing ultra cold atoms. The MOT is a vapor-loaded 
trap formed in a vacuum chamber with a base pressure 
~ 10~^ Torr. The six MOT beams are derived from an 
external cavity diode laser (ECDL) with the grating ar- 
ranged in Littrow configuration. The master diode laser 
is frequency locked to a saturation absorption feature 
produced in a Rb vapor cell. The laser power is increased 
by injecting the output into a slave laser, thus providing 
more than 20 mW of light in trapping laser beams of cross 
sectional area ~ 2 cm^. The slave output is switched 
and spectrally shifted as required with an acousto op- 
tical modulator (AOM) to a frequency about 18 MHz 
below the F — 2 F' ~ 3 trapping transition. 

The repumper laser is also an ECDL of the same design 
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as the MOT laser, and is locked to the = 1 ^ i^' = 2 
hyperfine transition. The repumper delivers a beam of 
maximum intensity ^0.6 mW/ cm?' and is delivered along 
the same optical path as the trapping laser beams. Re- 
pumper switching is controlled with an AOM. 

The cold atom sample is initially produced in the 
higher energy F = 2 hyperfine level. Direct absorption 
imaging measurements of the peak optical depth on the 
F ~ 2 ^ F' — 'i transition yielded, for this sample, 
~ 10 in a Gaussian radius of = 0.45 mm [ilJl. This 
corresponds to a total number of atoms ~ 3.2 • 10^ and 
a peak number density ^ 2.2 ■ VS^'^ atoms j err? . How- 
ever, the main goal is to transfer the trapped atoms to 
a carbon-dioxide (CO2) laser-based optical dipole trap. 
The 100 W CO2 laser operates at 10.6 [im. The laser is 
focussed to a radial spot size of ^ 55 /itm, and a corre- 
sponding Rayleigh range of ~ 750 /xm. The CO2 laser 
focal zone is overlapped with the MOT trapping region, 
while application of the CO2 beam itself is controlled 
by a 40 MHz AOM. After the atom sample is formed in 
the MOT, the CO2 laser is applied and the sample is 
compressed and loaded into the quasi static dipole trap 
(QUEST). This is done by detuning the MOT master 
laser 60 MHz to the low frequency side of the trapping 
transition, while simultaneously lowering the repumper 
intensity by an order of magnitude. The resulting tem- 
poral dark spot MOT loads the atoms into the lower en- 
ergy F = 1 hyperfine component, with about 15 % of the 
MOT atoms transferred to the QUEST. It is important 
to note that this transfer efficiency is measured after a 
QUEST holding period of about 1 second, during which 
the atomic sample collisionally evolves towards thermal 
equilibrium. Measurements of the QUEST characteris- 
tics, after the hold period, by absorption imaging, para- 
metric resonance, and the measured number of atoms 
transferred show a sample with peak density about 6 • 
10^^ atoms j cnt' and a temperature of 65 /iif. The 1/e 
lifetime of the confined atoms is greater than 5 s, lim- 
ited by background gas collisions. The residual magnetic 
field in the sample area, when the MOT quadrupole field 
is switched off, is estimated to be a few mG. 

Here we are concerned with parametric excitation of 
atoms confined to the QUEST and the sample character- 
istics. The sample is excited by modulation of the CO2 
power which, for a Gaussian focused beam, determines 
both axial and radial harmonic trap frequencies. Here 
these are about 1.25 kHz for radial excitation and 50 Hz 
for axial excitation; imposition of acoustic modulation on 
the CO2 laser AOM in the range 0-10 kHz is then suf- 
ficient to drive the fundamental parametric resonances. 
The characteristics of the modulation are the modulation 
frequency f, the modulation time T, and the modulation 
depth h. Detection of the result of excitation is made in 
two ways, each based on measurements on an image of 
the QUEST following a 3 ms period of free expansion. 
The expansion period allows the density of the sample 



to be reduced sufficiently that the QUEST becomes op- 
tically thin, so that measurement of light scattered from 
the sample is proportional to the number of atoms in the 
sample. In the first more traditional method, loss of a 
certain number of atoms is made through measurement 
of the total intensity of light scattered from the sample. 
This trap loss method measures the survival probability 
of atoms in the trap. In the second method, and the one 
we focus on here, the average peak intensity in the central 
region of the image is measured; this is proportional to 
the survival probability of atoms more localized spatially 
in the harmonic region of the trap, even after a time T of 
parametric resonance. The main measured quantities are 
fiuorescence images of the expanded atomic cloud, these 
being recorded for different h, T, and f. Characteristic 
results are shown in Fig. 2, where the images show a 
clear increase in the axial and radial Gaussian radii and 
a loss in peak intensity as T is increased at fixed h and f = 
2.5 KHz. We also draw attention to the change in shape 
of the cloud with increasing T. This is due to heating 
of the atom sample while it is confined, and consequent 
expansion in the weakly confining axial direction. 
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FIG. 2: Typical images of the expanded cloud of ultracold 
Rh atoms, after a period of parametric excitation. In all 
cases the cloud is permitted to expand for 3 ms prior to imag- 
ing. The quantity T is the modulation time, while Jo is the 
peak intensity in the central region of the image. To is the 
Gaussian radius of the image in the radial (vertical) direc- 
tion, which provides a spatial scale for the images. 

We begin our analysis by showing in Fig. 3(a) the 
excitation spectrum in the spectral vicinity of the funda- 
mental radial resonance. Two sets of data are displayed. 
In the first, we show the total number of trap atoms sur- 
viving the parametric excitation, while in the second we 
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FIG. 3: (a) Fundamental radial parametric resonances for h 
= 0.15 and T — 200 ms. (b) Fundamental axial parametric 
resonances for h = 0.20 and T — 200 ms. 



show the depletion of signal from the central spatial re- 
gion of the atom cloud image. The differences between 
the two are striking. First, even for the quite large mod- 
ulation depth of 15 % the total trap loss, as a fraction of 
the entire signal level, is clearly weaker than the depletion 
signal. The second main feature is that the trap loss sig- 
nal is shifted by -0.3(1) KHz to lower frequencies. Both 
aspects of the total trap loss signals represent significant 
experimental disadvantages if one wants to characterize 
the lower energy portion of the optical dipole trap poten- 
tial. We attribute the frequency shift to the anharmonic 
nature of the trapping potential; as the average amount 
of excitation energy of atoms in the trap increases, the 
resonance harmonic frequency for those atoms decreases, 
shifting the resonance to lower values. Because the res- 
onances can be intrinsically quite broad, this shifts the 
overall response to lower frequencies. Roughly speak- 
ing, one can envision the process as parametric evolution 
of an atom distribution localized deep in the trap to a 
warmer one which is more broadly distributed over the 



attractive part of the optical dipole potential. This leads 
to a 'freezing out' of the parametric resonance condition 
and to a depletion of the number of atoms in the deepest 
part of the atom trap. We emphasize that similar results 
are obtained for a range of h and T, but that the shift 
of the resonance to lower frequencies becomes more pro- 
nounced as h is increased. Within the range of our data, 
the resonance positions are insensitive to T. We illustrate 
the point further in Fig. 3(b), which shows the main ax- 
ial resonance for the same two measurement approaches 
used in Fig. 3(a), but for a larger modulation depth of 20 
%. In this data, there is no clear resonance for the trap 
loss signal, but a distinct one appearing at about 105 
Hz for the depletion signal. The spectral location of this 
resonance is, within the experimental uncertainty, where 
it is expected based on the spectral location of the fun- 
damental radial mode at 2.5 kHz, and assuming a trap 
formed at the focus of a Gaussian trapping beam. 

We now turn our attention to the role of the modula- 
tion time T, which represents the amount of time that the 
trap depth is modulated before the atoms are released. 
Physically, T represents the amount of time that energy 
may be delivered to the ultracold atom sample. We point 
out that this time scale is short compared to the 5 s hold 
time of the atom sample. The variation with T of the in- 
tensity at the peak of the image, for several modulation 
frequencies, is shown in Fig. 4(a). There we see that the 
peak image intensity appears to decrease exponentially 
for a range of f. This decrease is evident from the small- 
est T = 5 ms, and continues for a factor of 60 in T. In 
Fig. 4(b) the variation of the total trap loss with T for 
the same modulation frequencies is shown. In contrast to 
the peak intensity, the total intensity decreases approx- 
imately linearly with T, and further shows a threshold 
below which minimal variation is measured. We quali- 
tatively interpret these results as follows: the atoms in 
the trap are initially in thermal equilibrium at a temper- 
ature of approximately 65iJ,K. While the modulation is 
applied the average energy of the atoms is expected to 
increase. Further, the elastic collision rate is quite 
large at the experimental conditions, so the atomic cloud 
is expected to stay in approximate thermal equilibrium at 
some temperature. This temperature should increase as 
T increases. Even if no atoms were removed from the trap 
by this process, the peak intensity would still decrease, 
for both the spatial volume occupied by the atom cloud, 
and the rate of free expansion are increased as temper- 
ature is increased. However, as the trap depth is finite, 
some of the atoms have sufficient energy to leave the trap. 
This leads to a decrease in the integrated image inten- 
sity (due to trap loss), and a corresponding additional 
decrease in the peak image intensity. These arguments 
also explain why there is a threshold in the appearance 
of trap loss; the average amount of energy deposited per 
atom in the system by parametric excitation must be suf- 
ficiently large that a significant number of them have a 
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FIG. 4: (a) Peak survival probability as a function of T for 
several different parametric resonance excitation frequencies. 
Note the exponential decay with time for longer modulation 
times, (b) Relative survival probability as a function of T for 
several different parametric resonance excitation frequencies. 
Note the linear decay with time for larger T. h = 0.15. 



total energy in excess of the trap depth before trap loss 
is observed. Additional measurements of the Gaussian 
width of the images (as in Fig. 2) with increasing T pro- 
vide a direct measure of about 800 /iK/s for the heating 
rate, confirming this interpretation. We also point out 
that there is an important relationship between the re- 
sults of Fig. 4(a) and 4(b): when the decrease in the 
total number of atoms in the trap is accounted for by the 
data in Fig. 4(b), the remaining decrease seen in Fig. 
4(a) scales quantitatively as the mean cloud temperature 

3 /2 

as is expected for the temperature dependent de- 
crease in the peak intensity of a Gaussian atom cloud. 

Finally we point out that increasing T does not indefi- 
nitely increase the mean temperature of the atom cloud. 
As predicted in 13, above a certain critical amount 
of energy deposition in the system the temperature of 
the atom cloud remains constant, the additional energy 



going into removing additional atoms from the trap. As 
we show elsewhere [23|, this interpretation is supported 
by measurements on the variation of the Gaussian radius 
of the atomic cloud with T after ballistic expansion. The 
radius is seen to initially increase with T, but then to 
reach a steady value corresponding to a cloud tempera- 
ture of about 35% of the trap depth, consistent with the 
results of Gehm, et al. [l^]- 
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